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1 INTRODUCTION 



An activity of many compact objects - Active Galactic Nu- 
clei (AGN), Young Stellar Objects (YSO), microquasars - 
is associated with the highly coUimated jets. These jets are 
thought to be a natural outlet of an excess angul ar momen- 
tum of a central object and accreting matter (jHevvaertsI 
1 19961 ). T he latest observat ions indica ting the jet rotation 
in AGN (|Young et alll2007l ) and YSO (|Bacciotti et aj|2007l ) 
support this idea. The most attractive model for such out- 
flows is the MHD one dHeyvae rts 1996; Bland ford fc Pavnd 
1 19821 : iPelletier fc Pudridll992i ) . 

Of course, the main question w ithin this model 
is the coUimatio n itself felandfor d fc Pavnd 19821 : 
IPcUetier fc Pudrit j Il992l: ISautv fc fsinganosi 



ABSTRACT 

We demonstrate that the ID cyUndrical version of the Grad-Shafranov equation is 
more rich than classical self-similar ones, and more suitable for the astrophysical jets 
we observe. In particular, it allows us to describe the central (and, hence, the most 
energetic) part of the flow. Both relativistic and non-relativistic versions are discussed. 
It is shown that taking into account the finite pressure of the external media one can 
determine the magnetic flux within the central core. We found as well that for non- 
relativistic flows which are magnetically dominated near the origin the solution can 
be constructed only in the presence of the oblique shock near the base of a jet where 
the additional heating is to take place. 

zero electric current ||App1 fc Camenzind|[l993l ). but it can 
be constructed if the angular velocity vanishes at the jet 
boundary and if the exte rnal pressure is not equal to zero 
(Bcskin 1997; Beskin fc M alvshkin 200(1). 

Another result, obtained f or both relati vistic and non- 
relativistic cylindrica l flows (IChiueh et all 



1 19931 : lBogovalov|[l995l : iHevvaerts fc NormanI 



19911: lEichlerl 



20031 ). is that 



the poloidal magnetic field Bp has a jet-like form 
Bo 



B 



1 + rVrI 



(1) 



iShu c t al Il994l : lOuved fc Pudritj Il997l ). We 



1994 



, here that the coUimation is d ue to a finite exter n al gas 
and /or magnetic pressure (lAppl fc Came nzindl 1 19931 : 
iLerv et al|[l999l : iBeskin fc Malvshkinl l2000). Indeed, propos- 
ing that it is the external magnetic field Bcxt ~ 10~® G 
that plays the main role in the coUimation, we obtain 
''jot ~ Rin (Bin/Bcit)^^^ ■ Here r is the distance from the 
rotational axis, and the subscripts 'in' correspond to the 
values taken in the vicinity of the central object. The 
similar evaluation can be obtained for the external pressure 
Pext ~ B^^JStt. As for YSO Bin ~ 10^ G and Rin ~ i?©, we 
obtain rjct ~ 10^^ cm, in agreement with the observational 
data. Accordingly, for AGN {Bin ~ 10* G, Kin ~ 10^=* cm) 
we have rjet ~ 1 pc. It means that the external media may 
indeed play an important role in the coUimation process. 

The internal structure of cyl indrical jets was consid- 
ered b oth for non-relativist ic ( Gontopoulos fc Lovelacel 

[Hii; 



where rcorc ~ i'in7in/^- But this relation corresponds to a 
very slow (logarithmic) growth of the magnetic fiux func- 
tion: ^'(r) oc Inr. It means that if the jet core contains only 
a small part of the total magnetic flux the jet boundary 
is located exponentially far from the axis, with the mag- 
netic field being too weak to be in the equilibrium with the 
external pressure. In what follows we try to resolve this con- 
tradiction. 

Thus, we consider the following model: the fiow crosses 
all the critical surfaces while the effects of the external media 
are negligible. It allows us to use standard values of integrals 
of motion. As the supersonic wind expands, its pressure be- 
comes comparable with the external gas and/or magnetic 
pressure. The interaction of a fiow with external media re- 
sults in well coUimated jet which can be described by ID 
cylindrical equations. 



^Heyvacrts fc Norman 200^) and relativistic 

llChiueh et al 1991; Appl fc Camenzindl Il993l: lEichlerl 



2.1 



BASIC EQUATIONS 
Relativistic flow 



19931: lBogovaiov|[l9 96: Isto min fc Parievlll996l : lBeskinlll997l : 



Beskin fc Nokhrinai2006i ) fiows. In particular, it was shown 



For a cylindrical fiow one can write down electric E and 
magnetic B fields as well as the four- velocity of a plasma u 
in the standard form 



that for the constant angular velocity of plasma D,f it 
is impossible to obtain a reasonable solution with total 



B 



27rr 



21 



E : 



2ttc 



(2) 
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u = -B +7 

n 



c 



(3) 



Here n is the concentration in the comoving reference frame, 
and 7^ = u'^ + 1 is the Lorentz-factor. In other words, it is 
convenient to represent all the values in terms of a magnetic 
flux and a total electric current I, the angular velocity of 
plasma fip and the particle to magnetic flux ratio rj being 
constant on the magnetic surfaces: flp — flp{^), rj = J7(*I'). 
Accordingly, the tra ns-field Grad-Shafrano v (GS) equation 
can be rewritten as (|Beskin fc Parievlll993l ) 



1 d 
r dr 



, r dr / c 



G 



d^F 

2 Vdry "di" ' r27W2c4d* 
dr; 
d* 



M? 'd* " d* °' 

Here the entropy s = s(\[') is the fifth integral of motion, 

G^r'^{E-Q.YLf ^JvCL^c -M^r^E^, (5) 

A = is the Alfvenic factor, A^^ = A-R^mf /n 

is the poloidal Alfvenic Mach number, /i = mp(? + rripVj is 
the relativistic enthalpy, and the derivative d/d^ acts on the 
integrals of motions only. Finally, the relativistic Bernoulli 



equation Up = 7 

647r4r2 \dr ) 
where 



K 



1 has a form 



2 2 



(6) 



(7) 



and e' — E ~ QpL. Both equations contain relativistic inte- 
grals of motion 

= jfiT^c^ + L(*) = ru^firjc + ^, (8) 

which, as all other invariants, are to be determined from 
boundary and critical conditions. E.g., for the inner part of 
a fiow VE" ^ ^0 with a zero temperature one can choose 
nF(*) = fio, = Vo, and 

= MWincV 0*, L(^) = -^<K. (9) 
Multiplying now equa tion ([4t by 2 Ad\&/dr and using 



equation one can find l|Beskinl 119971 ) 



(e')= 



- 1 



22 2 ■ 



+ 



c 



+M 



^^3^2 ^2 g2 

;' d* de' 



d.M^ 
dr 

{e'f 



Afi'^r/^c^ 



1 - 



d^ d^F 
dr d* 

„ ^ , d* 1 dr7 
+ 2^4 — i 



(10) 



n \ OS / , 



+ 1- 



^Fr' 



c / dr ?7 d"!/ 

■M^ dl' da 
fi dr d*' 



where T is the temperature and Cs <C c is the sound velocity. 
Together with the Bernoulli equation (O it forms the system 
of two ordinary differential equations for the Mach number 
M^{r) and the magnetic flux ^'(r) describing cylindrical 
relativistic jet. Clear boundary conditions are 



*(0) = 0, P(rj,t) 



(11) 



where P — /Siv +p is the total pressure. Determining the 
functions A4^{r) and *I'(r), one can find the jet radius rjct as 
well as the profile of the current /, the particle energy, and 
the toroidal component of the four-velocity using standard 
expressions 

/ _ L~ npr^E/c^ 
2^ ^ ~ 



7 = 



1 {E -QfL) - M'^E 



(12) 



1 {E - Q.FL)nFr'^/c^ - LM'^ 



jj.r]rc 



1 - f7|rVc2 - 



(13) 



(14) 



2.2 Non-relativistic flow 

In the non-relativistic limit the electric and magnetic fields 
are determined by general expressions ((2]). On the other 
hand, equation (|3]) can be rewritten as 



Vn 



B + f^Fre^ 



(15) 



where pm = rripn is the mass density, and r?n(^) is the non- 
relativistic particle to magnetic fiux ratio. Accordingly, non- 
relativistic fiuxes of energy _En and z component of the an- 
gular momentum Ln are 

(16) 



(17) 



Further, algebraic relations (|12p . (|14|l can be rewritten as 



2^ = 



1-7W2 ' 

1 Q.ft'^ - -LnX^ 

= 

r 

where now 

47177,2 



.M2 



Pm 



(18) 
(19) 

(20) 



As a result, the non-relativistic Bernoulli equation 



(fiFr^ - L^M^f 



- 2r^nF 



(1-^2)2 1-^2 ' 

together with the non-relativistic limit of equation (flO 

dX^ 



(21) 



2e„ - 2it) nlr^ - (1 - VW^)cs 



dr 



>|2(2X2 - 1) 
+M"^^ + M'^r'^Q.v- 



l-M'^r^ 1-7W2 
^2d* de; 
dr d* 



d* dST^F 



+2 



w + 



(l-A^')cf 



dr d* 
A^d* dr?n 
77n dr d*? 



(i-x^)^(f:) +^ 

Pm \OS J TTlp 



d* ds 
dr d*' 



(22) 



where e-n = En — QpL^, determine the structure of a non- 
relativistic cylindrical fiow. 
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3 ADVANTAGES 

Certainly, the approach under consideration is one- 
dimensionaL For this reason, it has some properties similar 
to the other self-similar ones. In particular, one can eas- 
ily check that the singularity on the fast magnetosonic sur- 
face is absent. On the other hand, singularity appears on 
the cusp surface where the factors in front of dA4^/dr in 
UlUp and (|22p vanish. Nevertheless, in our opinion, this one- 
dimensional approach has some clear advanta ges in compar- 
ison with the standard self - similar ones ([Blandford fc Pavnd 
19821: iPelletier fc Pudrit3 Il992l : ISautv fc Tsinganod 1 1991 : 
Shu et al|[l99i ). 

First of all, it allows us to use any form of the five inte- 
grals of motion. Indeed, the self-similarity of a flow demands 
deflnite dependence of invariants which may not correspond 
to the real boundary conditions. E.g., for relativistic self- 
similar flo w the a ngular velocity Sip js to have the form 
f2p oc (iLi. Chin ch & Bcgolman 1992). It does not corre- 
spond neither to the homogeneous stellar rotation, nor to the 
Keplerian disk rotation. Moreover, this dependence has the 
singularity at the rotational axis. Thus, the standard self- 
similar approach cannot describe the central (and, hence, 
the most energetic) part of the flow. 

Further, classical self-similar approach cannot describe 
the region of the electric current closure. Finally, for the rel- 
ativistic magnetically dominated flow it is more convenient 
to use first-order equation (|10|) instead of second order GS 
equation for which it is necessary to be careful in taking into 
account small but important terms ^ 7~^. Indeed, the force 
balance equation (|10|l does not contain the leading terms 
PoE and j x B/c as they are analytically removed using 
Bernoulli equation. As a result, as 



\peE+j X B/c| 
|j X B/c| 



(23) 



4 INTERNAL STRUCTURE OF 
CYLINDRICAL JETS 

4.1 Relativistic cold flow 

4-. 1.1 General properties 

The solution of equations ([6| and (|10p for relativistic cold 
flow Cs = 0, s = depends essentially on the Mach number 
on the rotational axis Ml = M^(0) jBeskin fc MalvshkinI 
l200d ). For Ml > Ml,, where 

Mi = Til 
we have 



(26) 



M'' = Ml{ 1 



(27) 



the poloidal magnetic field corresponding to jet-like solution 
(O. Here i?L = c/S1f(0). On the other hand, for Ml < Ml, 
we obtain 



2Mla \RlJ 



M' = 



Ml 



Here 



(28) 



(29) 



(30) 



is the Michel magnetization parameter (jMichelll 19691 ') (7 a; cr 
for particle dominated fiow M^part ~ Wem)- It means that 
Bp ~ const, i.e., the solution has no jet-like form. 

As was already stressed, the solution p7p cannot be 
realized in the presence of the external media. Hence, one 
can conclude that for any finite external pressure Pext mag- 
netic field in the center of cylindrical jet Bo = A-Kr/fi'yin/Ml 
cannot be much smaller than Bmin = 47rj7^7in/A^cr- It gives 



all the terms in equation (|10[) are of the same order. 

In particular, in the limit r S> ^core, M^ ^ 
1 equation (I10|) can be rewritten in the simple form 
l|Beskin fc Malvshkinll2000l) 



d^ [ M^ 



M^ 



HrjflFr^iQlr^/c'^ + M^)' 



(24) 



Without the last term cx L^(\[') equation (|24p results in the 
conservation of the value H 



H 



^Iprjr 
M^ 



= const, 



(25) 



was found by iHevvaerts fc NormanI l|l989t ) for the conical 
magnetic field. It is the conservation of H that results in the 
jet-like solution ([ij. Indeed, as »7(^') « const and f2F(^') ~ 
const in the very center of a jet, we obtain M'^ oc r^. Using 
now the definitions M^ = A-Rff^/n and nup = r]Bp (and 
the condition Up ~ const in the very center of a flow), we 
return to ([l]). But, as we will see, the term containing 

i which appears to be missed previously) can be important 
Beskinl ll997). It is this term that can change the jet-like 
structure in a relativistic case. 



Bmin= B(i?L), (31) 

0-7in 

where B(J?l) = ^o/tt-Rl- 



4- 1-2 Central core 

Thus, for the external magnetic field Bext > -Bmin the in- 
ternal structure of a relativistic jet is to be described by 
relations (|28|l - (|29p when Bp ~ Bcxt. On the other hand, 
for Bext < Bmin iu the center of a flow (i.e., for r < 
7in-RL) the core with Bp « Bmin is formed. As was found 
(Beskinfc Nokhrinalliooel ). for a-^B{Ri^) < Bext < Bmin 
(and for r ^ ymRh) the solution can be presented as 

^ocr", M^ocr'', (32) 

the sum being a + b — 3. E.g., for Bcxt = Bmin we have 
a = 2, & = 1 (cf. (IIHI)-®), and for B^xt = g''^B{Ri^) we 
have a = 1, 6 = 2. 

The results presented above were reproduced recently 
both analytically and numerically. In l|Beskin fc Nokhrinal 
i2006i) it was shown that ID approximation becomes true for 
the paraboloidal outflow at large distances from the equa- 
torial plane z ^ g^^'^Rij where the flow becomes actually 
cylindrical. Up to the distance z — a^inRh one can use the 
relations (|28p - (|29|) . so the poloidal magnetic field does not 
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depend on r. The region z > aji^Rh corresponds to core- 
like solution (|32|l . Nevertheless, the transverse dimension of 
a jet remains parabolic: r\et <x z^^^ . Numerically the scalings 
H32p were confirmed bv lKomissarov et 

Remember that the existence of a cylindrical core 
with Tcore ~ 7inflL was predicted in many papers 
l|Hevvaerts fc Norman|[l989l : lBogovalovlll996r i. but the mag- 
netic flux 



mm 



(33) 



inside the core was unknown up to now. As we see, in the 
relativistic case the central core contains only a small part 
of the magnetic flux: 



Tin 

a 



(34) 



Nevertheless, as 6 > 0, such core-like flow can exist in the 
presence of the external media. 



4.1.3 Bulk acceleration 

As on the fast magnetosonic surface (rp ~ a^^^Rh) 
in the region of the diverging magnetic fi eld lines the 
bulk plasma Lorentz-factor yir-p) = cr^^^ l|Michell 1 19691 : 
iBeskin. Kuznetsova fc Rafikovlll998h . and, hence, here 



part 



^ a-'/' « 1, 



(35) 



an additional particle acceleration is possible as the trans- 
verse dimension of the diverging flow becomes larger than 
rp. Using equation (|13|) and the relation a + b = 3 one can 
flnd that in the whole region _Bext > B{Ri^) (z < a^Rh 
for the parabolic flow) the Lorentz-factor can be determined 



as 

7 « r/Rh. 

Accordingly, one can write down 

1/2 



VKpart 



B{Rl) 



B,: 



(36) 



(37) 



It means that for _Bext ~ o" ^B{Ri^) {z ~ (J^Rl for the 
parabolic flow), where the transverse jet dimension 



(38) 



almost the full energy transformation from the Poynting flux 
to the particle energy flux can be realized. In particular, for 
the particle moving along the parabolic magnetic fleld line 
one can obtain 

jiz) ^ {zlR^fl\ (39) 

This scaling was confirm ed numerically as well (jMcKinnevI 
l2006l : lNaravan et alllioo^ ). 

It is necessary to stress that relation (|36|l takes place 
only if one can neglect the curvature of the magnetic sur- 
faces. Indeed, for the magnetically dominated case in the 
limit r 2> rp the leading terms in two-dimensional GS equa- 
tion PeE -I- J X B/c can be rewritten in the simple form 
jBeskin fc Nokhrina|[2006l ) 



in.V(B^)-^(n.R. 



Bi - 



(n • er) = 0. 



(40) 



Here R^ is the (poloidal) curvature radius of magnetic sur- 
faces, Rc is the unit vector in the direction of curvature 
radius growth, and n = V'I'/|V*I'|. 

Neglecting now the curvature term and using standard 



relations B^ ~ Bpv/Ri^ and B^ 



E^ 



B^/j resulting 



from ((5)) and we return to H36p . On the other hand, if 
the curvature is important, then one can neglect the first 
term in (|40)) . and we obtain (jBeskin. Zakamska fc Solll20o3 ) 



7 « (Rc/r) 



1/2 



(41) 



This scaling taking place for the split-monopole geome- 
try outside the fast magnetosonic surface correspond s to 
(|Tomimatsulll994l : [Beskin. Kuznetsova fc Rafikovll 19981 ) 

7 « a'''='ln'''='(r/rF). (42) 

Remember that for r < rp we have a "linear" acceleration 



As was demonstrated numerically (|Naravan et al|[2006l ) , 
it is the parabolic fiow that terminates these two asymptotic 
solutions. If the magnetic surfaces have a form z oc then 
for k > 2 (when the coUimation is even stronger than for a 
parabolic fiow) one can use the relation 7 « r/Ri^. On the 
other hand, for 1 < fc < 2 at large distances the particles 
acceleration is not so effective, so that 7 « [Rc/rY^'^. As 
Rc ~ {z )^ jz" for z' ^ 1, where z' = dz/dr, we obtain for 
the Lorentz-factor of a particle moving along the magnetic 
field line z = zir) 

7 oc zf'^-i)''*, (43) 

in full agreement with the cases fc = 2 (|4ip and k = 2 (|42l) 
considered above. Accordingly, the total energy transforma- 
tion can be realized if the jet width is 

roq ~ a^/C'-^'i^L. (44) 

For fc = 3/2 these sca l ings w ere confirmed numerically by 
iBarkov fc Komissarovl (120081 ). For fc > 2 the evaluation 
Teq ~ cr _Rl is to be used. Thus, effective particle acceler- 
ation can take place only if rjet ^ (jRh , and if the curvature 
of magnetic surfaces is not important. 



4.1.4 the center of the self-similar solution 

The approach under consideration allows us matching the 
self-similar solution to the rotational axis. Indeed, for 
the self-similar invariants |Li. Chiueh fc BegelmanI Il992l : 



IContopoulos fc Lovelacelll994 ) 

np(*) = f7o(*/*b)"'', 

£(*) = £o(*/*b)'"'"' 

L(*) = Lo(*/*b)'"'', 
VW = ^o(*/*b)'-"', 



(45) 
(46) 
(47) 
(48) 



one can seek the solution of two-dimensional GS equation 
for * > *b in the form *(_R, 9) = J?^/''e(e'), where R is the 
spherical radius. Hence, for 9 <^ 1 one can write down 



^(J?,^) =AR^^^9' 



(49) 



where A — const, so the cylindrical radius of the boundary 
^' — ^'b can be written as 



rb(z) 



i/a,T,l/o. 1-1/ab 



(50) 
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Figure 1. Internal structure of magnetically dominated relativis- 
tic jet. The short dashes line is an exponent a of the flux function 
'3/ I I49I I. The long dashes line represents the logA^^, and the solid 
lines represent an exponent bof H52| l for different values of the 
parameter /3. The thin dashes lines are the analytical exponents 
6 = 3 — 6/3. The non-dimensional radius x = Q,-p(0)'r / c. 



On the other hand, as was already stressed, in the cen- 
tral part of a flow ^I* < the self-similar approach cannot 
be used. For simplicity we assume that Q.y = f2o = const 
and ri — rio = const for 4' < "tb- Then far from the equa- 
torial plane where z S> r ^ 1) one can integrate ID 
cylindrical equations (|5J and pup considering z m R as a 
parameter. As a result, using solutions (|28|l and (|29[) . we 
obtain for M'^{z) = Al'^(rb) two different expressions for 
particle and magnetically dominated flows. 

For magnetically dominated How (rt S> ^inRh) we have 



Ml{z) = 



47r rjofi 



3-3/a0 



1-3/a 



(51) 



Besides, for 'i' > vpb (r > rb) one can seek the solution in a 
form 



As a result, equations (|6]), (|10p give 
a = 2, 6 = 3-6/3, 



(52) 



(53) 



the first relation demonstrating that poloidal magnetic field 
is to be homogeneous: Bp ~ const. Substituting now r = R9, 
we see that for 'I' > 'I'b 



M=ce, 

where the coefficient 

C oc Mb^{R)R''/rl{R) oc i?3+(6-3)/2/3^ 



(54) 



(55) 



in agreement with the self-similar property, does not depend 
on R. 

For particle dominated How (rt <C 7ini?L) we find 
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Figure 2. The same for particle dominated relativistic jet. The 
thin dashes lines are the theoretical exponents 6 = 2 — 4/3 



Ml 



_42/a^l-2/« 



(56) 



Seeking again the solution for 4* > ^tb in a form (|52|l . we 
have 



4/3. 



(57) 



This solution ensures independence of a coefficient C on R 
as well. The results of numerical integration of the system 
([6|, (jlOp for the cold flow are presented in Figs. 1, 2. As we 
see, there is very good agreement between numerical results 
an analytical asymptotic behaviour (|53|) and (|57|l . 



4.2 Non-relativistic flow 

4-2.1 Central core 

For the non-relativistic cold flow in the central part of a jet 
one can use general expressions @ 



Ln(*)=lio- 



(58) 



where Vin can be considered as a constant, and the non- 
dimensional longitudinal current io = j/jc} depends now 
on the angular velocity Qp. Here jcj = O.oB/2-r is the 
Goldreich- Julian current density. For fip ^ ficr, where 

1/2 



«in_ 
^in 



Pintail 

Bl/9, 



(59) 



corresponding to a particle dominated outfl ow near the star , 
the 2D problem can be so lved analytically l|Bogovalovlll995l : 
iBeskin fc Qkamotd [2OO0I 1 . and we obtain io = c/vi-n. For a 
magnetically do minated flow ne ar the origin one can write 
down (see, e.g., l|Lerv et a]||l999l ')') 
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-2/3 



10 



Win Vricr/ 

Nevertheless, for 'I' < ^i^i, where 
47r^wfn77o 



*in = 



(60) 



(61) 



the flow remains particle dominated: ~ ■t'hi/2. The con- 
dition = x]>o just corresponds to f2p = f2cr. Remember 
that in the non-relativistic case the flow can pas s smoothly 
the c ritical surfaces only if Wpart(f"F) ~ M4m(?'F) (|Heyvaertsl 
Il996l ). Thus, the flow at large distances is to be particle dom- 
inated. 

At first, let us consider the simplest case of sub Alfvenic 
cylindrical flow < 1. Solving equations (|21|l and (|22|) one 
can flnd that the poloidal magnetic field remains constant 
inside the jet up to the very boundary. Thus, one can put 
B(0) = -Bext. There is a simple physical explanation why 
the homogeneous poloidal magnetic field is a solution of the 
trans-field equation for a subsonic flow. The point is that 
for < 1 the energy density of the poloidal magnetic 
fleld Bp/Syr is much larger than both the energy density of 
the toroidal magnetic field /8tv and the energy density of 
particles puiV^/2. As a result, the trans-field equation can 
be rewritten as 



d 
dr 



Stt 



0. 



(62) 



Hence, for sub Alfvenic fiows the homogeneous poloidal mag- 
netic field is a solution of the trans-field equation for arbi- 
trary invariants En{^) and Z/n(^'). But such sub Alfvenic 
flow can exist only in the presence of large enough external 
magnetic fleld Bext > -B(rF), where 



(63) 



For ordinary YSO B{rF) ~ 10^ G, so sub Alfvenic jets 
cannot be realized. 

On the other hand, cylindrical trans- Alfvenic flow can 
not also be realized both for the center part of a flow or for 
the self-similar region. To proof this proposition we must 
make two assumptions. We suppose that the derivative of 
the Alfvenic Mach number remains finite at the Alfven sur- 
face (AS), i.e. L — fipi'^l^g = 0. We also assume that the 
total current / is not closing at the AS strictly. 

Let us suppose that the flow in the center of a cylindrical 
jet is sub- Alfvenic and is about to cross the AS: = 1 — e. 



5, where e > 0, 5 > 0. In this case one can 



write down the leading terms of equation (|22|l as 



(2en -I- nlr''' 



dr 



o2 2 



d* , ^ 



den 
d* 



2 d* 



+ 2e, 



1 dri 
T] d* 



(64) 



The term in r.h.s. part of equation is equal to zero for the 
inner part of the flow, or it is negative for the self-similar 
integrals. As the total current does not close at the AS, 



const 7^ 0, 



(65) 



i.e., S — 0(e), we can neglect the first term in r.h.s. bracket 
in 1(64} . This leads us to a conclusion that the Mach deriva- 



tive in the vicinity of the AS is negative. However, if we 
assume that should reach the unity, there must be at 
least one point in the vicinity of the AS, where the deriva- 
tive is positive. We have come to a contradiction, so the 
transition of the AS is impossible. One can easily prove by 
analogy that the trans- Alfvenic flow is impossible also if the 
flow is super-Alfvenic in the center. 

Thus, super Alfvenic cold cylindrical flow must be su- 
personic at the rotational a xis: Mq > 1 . In this case we 
return to the jet-hke solution (|Eichleilll993l : lBogovaTovlll995l ') 



(66) 



But, as it was already stressed, in the presence of a flnite 
external pressure this solution is possible only if the central 
core r < rcoro = Vin/^ contains almost all magnetic flux 
^'o. This can be realized only for a slow rotation Q,y <C f^cr. 
In this case the magnetic field on the axis cannot be much 
smaller than 



-^min 



(67) 



Integrating now equations (121(1 and (|22p . one can obtain that 

^min 



ln(l + B^in/B 

c 

Accordingly, 



ln(l + Bmin/B cxt ) 



(68) 



(69) 



This structure was reproduced numerically as well 
(|Lerv et allll999l ). 

But for fast rotation Qp ^ Qct the reasonable solution 
cannot be realized as the core magnetic fiux ^Pcoro is much 
smaller even than the fiux 'I'ln (|61|) within the central part 
of a flow. Indeed, as according to definitions p5|l and (|20|l 
one can write down 



Sp(0) - 

we obtain for TMo > 1 



to "in 

2cA/2 



< 



-2/3 



< 1. 



(70) 



(71) 



It means that the cold cylindrical fiow resulting from the 
interaction of supersonic, fast rotating wind with external 
media cannot be realized. 



4-2.2 Heating at the oblique shock 

To resolve this contradiction, one can propose that in 
the observed non-relativistic supersonic jets an important 
role in the force balance may play the finite temperature. 
E.g., additional heating can be connected with the oblique 
shock near the base of a jet (|Bogovalov fc Tsinganosll2005l : 
iBromberg fc LevinsonI [2007h . It is well known that such a 
shock is needed jto explain the emission lines observed in 
jets from YSO (jSchwartz 1983). This situation is alike the 
pure hydrodynamical supersonic outfiow meeting the wall. 
The hydrodynamical analogy is all the more reasonable as 
the non-relativistic supersonic outfiow is to be particle dom- 
inated. 
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To evaluate the thermal terms in equations (|2ip and 
(|22p we consider pure hydrodynamical shock wave turning 
the spheric ally symmetric sup ersonic flow into the cylindri- 
cal jet (see (|Beskin et alll200"9r ) for more detail). For the field 
lines of the pre-shock flow at the inclination angles to the 
rotational axis less than the critical one, we seek the shock 
position so as to turn the flow into a cylinder. The critical 
angle for the pre-shock sound Mach number i S> 1 is 
equal to 

6'n.ax = sin-'(l/7), (72) 

where 7 is the polytropic index. In particular, for 7 — 1.2 
we have Smax ~ 56°. For the rest field lines we model the 
shock position to transit smoothly from the 9 = Smax to the 
equatorial field line 6 = n/2. Knowing now the shock posi- 
tion, we can calculate the entropy jump As for every field 
line. The rest integrals of motion, according to conservation 
laws, are to be unbroken on the oblique shock. 

Thus, we can solve one-dimensional equations (|21|l and 
H22p taking into account the effects of the heating on a shock 
through the corresponding thermal terms. We find that for 
the super Alfvenic flow the scaling ^' oc In r holds no more, 
so the jet boundary is located at the finite distance from the 
rotational axis. Obtained jet parameters for typical TTauri 
star (a jet radius Rjct ~ 10^^ cm, a temperature behind 
a shock needed to give rise to the forbidden emission lines 
T ~ 10* K, a poloidal velocity Vp ~ 10^ -^3 • 10^ cm/s, and a 
toroidal velocity ~ lO" cm/s) are in agreement with the 
observational data. 
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Figure 3. Internal structure of particle dominated non- 
relativistic jet. The short dashes line is an exponent a of the flux 
function Nt. The long dashes line represents the logA^^, and the 
solid line represents an exponent b of for different values of 
the parameter /3. The thin dashes lines are the theoretical expo- 
nents 6 = 3 — 6/3. The non-dimensional radius x = Qp{0)r/vin. 



4-2.3 In the center of the self-similar solution 

For non-relativistic cold outfl ow the self-similar invariants 
are (|Blandford fc Pavne|[l983 ) 



nF(*) = no(*/*b)""'/', 

£„(*) = £o(*/*b)"'', 

L„(*) = Lo(*/*b)'''", 

77n(*) = r;o(*/*b)'~'''/'. 



(73) 
(74) 
(75) 
(76) 



Again, one can seek the solution of the two-dimensional GS 
equation for * > *b in the form 'i'{R,9) = _R^'''^e(6l), so 
that for e < 1 



*(i?, 



AR 



l/0oa 



(77) 



Then, the cylindrical radius of the boundary ^ = can 
be written as 



rb(^) = 



(78) 



As a result, integrating one-dimensional cylindrical 
equations (|2ip and (|22l) in the region ^' < 5'b for parti- 
cle dominated flow, i.e., using the solutions (|28p - (|29p , we 
obtain for Sip — Qq and 77 = 770 



(79) 



On the other hand, for 'i' > ^'b {r > rb) one can seek again 
the solution in a form A4^{r) = Mh^{r/rb)''- As a result, 
equations ([JU, (HU give 



4/3. 



(80) 



Substituting now r = R6, we see that 

X^=Cr. (81) 

Again, the coeflicient C oc Mb^{R)R'' /r^{R), in agreement 
with self-similar property, does not depend on R. 

The results of numerical integration of the system (|2ip , 
(|22p for the cold non-relativistic particle dominated flow 
are presented in Figs. 3. Here we also see very good agree- 
ment between numerical results and analytical asymptotic 
behaviour (|80p . The rapid growth of the exponent b close to 
the end of the calculation is due to drop of Mach number 
close to unity where the flow is close to the Alfvenic surface. 
However, as we have showed, the smooth transition of the 
AS is impossible for the chosen self-similar integrals, so we 
must stop our cal culation at this point. T he power 6 = 2 — 4/3 
was confirmed bv lMatsakos et all (|2008l ) as well. 



5 CONCLUSION 

Thus, the cylindrical Grad-Shafranov equation has definite 
advantages in comparison with the standard self-similar 
ones. Using this approach it was demonstrated that in the 
relativistic case an effective particle acceleration can take 
place only if rjct ^ (jRl, the curvature of magnetic surfaces 
playing no role. We found as well that for non-relativistic 
flows which are magnetically dominated near the origin the 
solution can be constructed only in the presence the oblique 
shock near the base of a jet where the additional heating 
is to take place. In all cases the magnetic flux within the 
central core was determined. As was demonstrated, for rel- 
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ativistic flow the central core is to contain only a small part 
of the total magnetic flux. For the non-relativistic outflow 
the situation is opposite. 
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